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Applying  variation-correlation  method  to  images  of  fractography  obtained  from  the  scanning  electron 
microscope  (SEM),  it  has  been  found  that  there  exists  a fractal  characteristic  length  within  which  the 
fractured  surfaces  are  fractal.  Investigation  shows  that  the  fractal  characteristic  length  can  represent 
the  statistical  maximum  size  of  texture  of  the  SEM  image.  Multi-magnification  fractal  analysis  has 
shown  that  fractography  cannot  be  described  by  a single  fractal  dimension  but  rather  a series  of 
fractal  dimensions.  Fractal  study  on  the  fractured  surfaces  of  gray  iron  with  different  grain  size  has 
shown  that  there  exists  a positive  relationship  between  tensile  strength  and  fractal  dimension. 
However,  no  essential  relationship  between  impact  toughness  and  fractal  dimension  for  HP26Cr35Ni 
alloy  could  be  obtained.  Fractal  dimension  is  sensitive  to  changes  of  mechanical  properties  caused  by 
geometric  factors  such  as  grain  size  and  is  suitable  to  quantitatively  describe  the  irregularity  of 
fractography.  For  a given  fracture  mode,  fractal  dimension  is  not  universal. 


1 Introduction 

Fracture  phenomenon  is  one  of  the  most  intriguing  areas  in  material  science  and 
engineering  because  it  is  associated  with  the  life  and  safety  of  engineering  parts.  There  is 
a lot  of  information  related  to  fracture  mechanism  in  fractography.  Fractal  based  studies 
on  fractography  become  important  because  of  the  problem  of  variance  of  conventional 
quantitative  parameters  with  scales.  During  the  last  twenty  years,  many  fractal-based 
methods  for  quantitative  analysis  on  fractography  have  been  developed.  The  methods  can 
be  categorized  into  the  following  three  classes:  (a)  slit  island  analysis  [1);  (b)  profile 
analysis  [1'2 * * 5];  and  (c)  3-dimensional  (3D)  surface  analysis  [6_8].  Among  them,  3D  surface 
analysis  is  the  most  interesting.  The  fractal  characteristics  of  fracture  surfaces  and  the 
relationship  between  fractal  parameters  and  mechanical  properties  have  been  investigated 
using  these  methods.  Some  results  justify  that  fractal  dimension  can  be  a measure  of  the 
toughness  of  a material [1, 9’ 10].  On  the  other  hand,  some  results  have  shown  the  existence 
of  a universal  value  of  fractal  dimension  for  a given  fracture  mode.  It  has  been 
conjectured  that  the  value  of  fractal  parameter  a (a=3-D  where  D is  fractal  dimension)  is 
0.8  for  ductile  fracture  [11]  and  0.87  for  brittle  fracture  [3l  These  values  imply  that  there  is 
no  correlation  between  fractal  dimension  and  toughness  because  a universal  value  denies 
any  dependence  of  fracture  features  on  a specific  toughness.  In  this  paper,  a new  3D 
surface  model  based  on  variation  profile  analysis  has  been  introduced.  Using  this  model, 
the  fractal  characteristics  of  different  fracture  modes  and  the  correlation  between  fractal 
dimension  and  mechanical  properties  have  been  studied. 


2 Experimental  Procedures 

A typical  ductile  fracture  of  microvoid  coalescence  was  obtained  by  tensile  testing  a 

specimen  of  Assab  760  medium  carbon  steel.  A transgranular  cleavage  mode  was 

obtained  by  fracturing  a Charpy  V-notch  specimen  of  the  same  steel  after  immersing  it  in 
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the  liquid-nitrogen  for  60  minutes.  Fracture  of  a pure  AI2O3  ceramic  specimen  appeared  to 
be  brittle  intergranular  mode.  To  analyze  the  influence  of  geometric  and  non-geometric 
factors  on  fractal  dimension  of  fractography,  specimen  A of  gray  iron  was  modified  with 
0.4%wt  75SiFe  inoculation,  while  specimen  B,  with  0.3%wt  FeSiCuMo  metal  stream 
inoculation  to  improve  the  graphite-eutectic  microstructures.  Tensile  strength  ot,  was 
obtained  by  tensile  testing  the  various  materials.  A series  of  specimens  of  HP26Cr35Ni 
alloy  with  the  size  4x4x35(mm)  were  taken  from  same  regions  of  solidification.  Thermal 
shocks  before  mechanical  testing  were  carried  out  0,  5,  10,  15  and  20  times.  In  the  thermal 
shock,  the  specimens  were  heated  to  1 100°C  for  10  minutes  followed  by  water-quenching. 
After  the  thermal  shock,  impact  toughness  was  obtained  by  impact  testing  without  Charpy 
V-notch.  The  nominal  compositions  of  these  materials  are  listed  in  Table  1.  All  fracture 
surfaces  were  imaged  by  SEM  followed  by  digitizing. 


Table  1.  Composition  of  materials  investigated  in  wt%. 


Materials 

C 

Si 

Mn 

P 

S 

Cr 

Ni 

Fe 

760  steel 

0.5 

0.25 

0.5 

Bal. 

3.0 

0.12 

0.7 

<0.3 

<0.05 

Bal. 

HP26Cr35Ni 

0.55 

2.5 

2.0 

<0.04 

<0.04 

26 

35 

Bal. 

3 Variation-Correlation  Method 

The  variation  method  has  been  introduced  by  Dubuc  et  al  (4].  In  this  method,  an  e- 
oscillation  for  the  height  function  h(x)  of  a profile,  v(x,  s),  is  defined  as  follows: 

v(x,f)  = /imax(x)-/jmin(x),  I * - x J<  <?,  (!) 

The  integral  V(e)  may  be  defined  as: 

V(s)  = jv(x,s)dx.  (2) 

V(e)  is  termed  the  s-variation.  The  limit  of  V(e)  is  zero  for  s— >0  due  to  continuity  of  the 
function  h(x).  The  rate  at  which  V(e)  scales  when  £—>0  is: 

F(£)CCE2-D,  (3) 

which  provides  a method  to  calculate  the  fractal  dimension  D.  The  variation  method  can 
be  interpreted  as  a numerical  technique  of  box  counting  [4].  It  has  been  demonstrated  to  be 
more  accurate  than  any  other  techniques  for  fractal  analysis  of  self-affine  curves  [4,  . 

Based  on  variation  method,  a variation-correlation  model  for  3D  surface  has  been 
developed  [8].  In  this  model,  it  is  assumed  that  H(x,y)  is  a dimensionless  field-like  variable 
of  a surface  which  may  denote  distribution  of  height  or  angle.  If  Hmax(x,y,  e)  and  Hmm(x.y,  £) 
are  the  local  maximum  and  minimum  of  H(x,y)  in  an  ^neighborhood  respectively,  a 
variation-correlation  function  VC0Xs)  is  defined  as  follows 


Kor  (e)  =<  [Hm„  (x,  y,  e)  - Hmm  (x,  y,  e)]2  > 


(4) 
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The  symbol  <...>  refers  to  two  consecutive  averaging  operations.  If  the  variation- 
correlation  function  Vcof£)  of  the  variable  H{x,y)  on  a surface  obeys  a power  law  in  some 
significant  range  of  scales  sm,„  < £ < emax,  i.e., 

rw(e)cce2“  (5) 

where  the  fractal  parameter  a = 3 -Dcor,  the  variation-correlation  is  called  fractal  variation- 
correlation  function. 


4 Results  and  Discussion 


4. 1 Fractal  Characterization  of  Different  Fracture  Modes 

The  typical  fracture  modes  of  microvoid  coalescence,  transgranular  cleavage  and  brittle 
intergranular  fracture  are  shown  in  Fig.l  (a),  (b)  and  (c).  Fig.l  (d),  which  gives  the  Vcor-s 
curves  of  these  fracture  surfaces,  shows  that  the  variation-correlation  appears  as  lateral 
growth  only  within  the  range  £ < £max.  When  e > the  Vcor-£  curves  appear  to  reach 
plateau  value.  These  observations  imply  that  the  following  short-  and  long-range 
relationships  of  variation-correlation  versus  scale  are  true 


£ « £ 


£»£„ 


(6) 


From  Eq.(6),  the  log-log  plot  of  Vcor  versus  £ should  be  a straight  line  in  the  short-range 
of  scales,  of  which  the  slope  is  different  from  that  in  the  long-range  of  scales.  To  estimate 
emax,  a linear  measure  Mp  is  introduced  as  follows 


Mf  = 


1 -5 
1 + 6 


(7) 


where  8 is  the  degree  of  disperse  expressing  the  fluctuation  of  the  set  of  points  (x,  y,)  from 
the  fit  line.  8 is  defined  as 


5 = 


I>,  ~x)(y,  -y) 

»= 1 

£(x,  -x)2 
1=1 


£(T/-t)2 

_M 

2>,-x)2 

i=i 


where 


■:  -7^5- '. 

;;  t '■  ’* 
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Microvoid  type 
Transgranular  type 
Intergranular  type 


Microvoid  type 
Transgranular  type 
Intergranular  type 


Fig.l.  SEM  fractography  and  fractal  characteristics,  (a)  Microvoid  fracture;  (b) 
Transgranular  cleavage  fracture;  (c)  Brittle  intergranular  fracture;  (d)  Variation- 
correlation  versus  scale  curve  based  on  variation-correlation  analysis;  (e)  Height- 
correlation  versus  scale  curve  base  on  height-height  correlation  analysis.  The 
resolution  of  the  images  is  68.75u.m/200pixels. 


129 


The  value  of  MF  falls  in  the  range  between  0 and  1 . When  MF  tends  to  1,  the  distribution 
of  the  points  (x,,  y,)  is  perfectly  linear.  When  MF  tends  to  0,  the  distribution  becomes  non- 
linear. Table  2 shows  the  values  of  D and  estimated  with  Mp= 0.7,  while  the 
correlation  coefficient  r of  log-log  plot  is  over  0.98,  indicating  Mp  is  a better  linear 
measure  than  r. 


Table  2.  The  estimation  of  fractal  parameters  and  the  maximum  size  dmm  of  the  basic 
micro-morphology  element  on  the  fractography 


Fractography 

Variation-correlation 

Height-height  correlation 

-^max 

^cor 

£max 

r 

A if 

A,.h 

^max 

r 

MF 

Microvoid 

type 

2.58 

64 

0.984 

0.7 

2.68 

5 

0.982 

0.7 

75 

Transgranular 

type 

2.54 

17 

0.983 

0.7 

2.58 

5 

0.972 

0.7 

31 

Intergranular 

type 

2.44 

26 

0.983 

0.7 

2.33 

10 

0.982 

0.7 

35 

The  above  results  confirm  that  there  exists  a characteristic  length  cmnx  so  that  the  log- 
log  plot  of  variation-correlation  is  linear  within  amax  for  the  fractography  of  a real  material. 
Different  from  an  ideal  self-afFine  profile,  for  example,  Weierstrass-Mandelbrot  curve 
whose  Vcor  appears  as  lateral  growth  at  all  possible  scales,  the  morphologies  of  fracture  of 
real  materials  are  constituted  by  those  basic  micro-morphology  elements  such  as  dimples, 
cleavage  facets  or  boundary  facets.  The  basic  micro-morphology  elements  are  affine  at  a 
limited  scale,  to  form  the  texture  of  SEM  images.  Assuming  that  Vmm  is  the  expected 
value  of  maximum  variation  of  the  texture  at  all  possible  scales,  Vmax  is  bounded  within  a 
scale  range  (s  < £max)  because  of  the  affinity  and  periodicity  of  the  fractography  textures. 
Therefore,  this  characteristic  length  is  associated  with  the  statistical  maximum  size  of  the 
fractography  textures.  When  the  measuring  scale  e is  greater  than  the  largest  and 
smallest  values  on  the  fractography  are  always  included,  resulting  in  Vcod,  £)=V2 max. 

From  the  data  shown  in  Table  2,  the  value  of  Dmr  for  ductile  fracture  is  the  largest, 
while  for  brittle  intergranular  fracture,  is  the  smallest.  The  results  are  well  in  agreement 
with  the  irregularity  of  these  fracture  surfaces.  Fig.l  shows  some  small  dimples  are 
located  inside  the  larger  ones.  The  fractured  surface  experiences  a large  amount  of  ductile 
deformation  and  hence  dimples  at  different  scales  are  formed.  Such  fracture  surface  has 
multi-scaling  similar  structures  leading  to  large  smax,  within  which  variation-correlation 
follows  a lateral  growth  with  the  increasing  scale.  It  must  be  noted  that  a large  emax  does 
not  necessarily  mean  multi-scaling  similar  structures,  because  it  is  a parameter  related  to 
the  maximum  size  of  the  surface  texture.  For  the  brittle  intergranular  fracture,  the 
fractography  is  constituted  by  smooth  boundary  facets,  resulting  in  a decrease  in  Dcor.  The 
irregularity  of  the  transgranular  cleavage  fracture  is  medium  in  spite  of  the  rather  rugged 
surface  because  of  the  smoothness  of  the  surface  of  each  grain. 

For  the  sake  of  comparison,  the  curves  of  height  correlation  versus  scale  using  height- 
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height  correlation  analysis  [1215]  are  also  shown  in  Fig.l  (e).  The  height-height  correlation 
function  C(e)  is  defined  as 


C(s)  =<  h(e)  - h(0)  > 

where  h( 0)  is  the  height  of  reference  point  and  h{s),  the  height  of  the  points  with  e 
distance  from  the  reference  point.  The  estimated  values  of  the  fractal  parameters  are 
shown  in  Table  2.  Fractal  dimensions  from  height-height  correlation  appear  the  same 
tendency  as  those  from  variation-correlation  in  characterizing  the  irregularity  of  the 
fracture  surfaces.  However,  the  characteristic  length  smax  in  the  height-height  correlation 
analysis  represents  the  correlation  length  !12'151  but  not  the  size  of  texture.  Consequently, 
the  fractal  characteristic  length  in  variation-correlation  analysis  is  a more  suitable 
parameter  in  quantitative  fractography  because  its  estimated  value  is  able  to  represent  the 
maximum  size  dm.ix  of  the  basic  micro-morphology  elements.  Consider  a simple  case  of 
brittle  intergranular  fractography,  which  simply  consists  of  uniform  morphology  elements 
with  maximum  diameter  dmax.  In  variation-correlation  analysis,  when  dm ax/2  < s < dmm, 
some  estimated  values  of  variation  in  statistical  measurement  have  already  reached  the 
value  of  Fmax.  When  s > dmm,  all  the  estimated  values  of  variation  are  Vmax.  Therefore, 
when  £ > dmsJ2,  the  linearity  of  the  log-log  plot  becomes  worse  so  that  the  estimated 
value  of  £max  falls  into  [dmJ2,  dmix\  while  its  fluctuation  depends  on  the  given  value  of 
A/f.  The  results  in  Table  2 hence  verify  the  validity  of  the  analysis.  Therefore,  the  fractal 
characteristic  length  can  be  employed  as  a statistical  measure  of  the  maximum  size  of 
texture  of  fracture  surfaces. 

4. 2 Relationship  between  fractal  dimension  and  mechanical  properties 

Fig.2  shows  the  SEM  fractography  and  the  multi-magnification  fractal  curves  for  the  gray 
iron.  Specimen  A has  been  modified  by  FeSiCuMo  metal  stream  inoculation.  Thus, 
microstructures  of  the  matrix  and  graphite  structures  were  improved,  leading  to  an 
increase  in  tensile  strength  <xb  (<rb=314MPa  for  sample  A and  255MPa  for  sample  B).  The 
elongation  for  gray  iron  is  almost  zero.  The  size  of  graphite  and  graphite-eutectic  cells  is  a 


Fig.2.  SEM  fractography  of  samples  A and  B of  gray  iron,  (a)  SEM  fractography  of 
sample  A;  (b)  SEM  fractography  of  sample  B.  The  resolution  of  the  images  is 
421pm/160pixels. 
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very  important  geometric  factor.  Small  grains  will  lead  to  multi-scaling  fractal  structures 
of  the  material  and  the  fracture  surface  becomes  more  complicated.  The  cleavage  facets  of 
specimen  A,  which  is  associated  with  grain  size,  are  smaller  than  those  of  specimen  B, 
leading  to  the  more  complicated  fractured  surface  as  shown  in  Fig.2  (a).  This  implies  that 
the  fractal  dimensions  of  specimen  A is  larger  than  that  of  specimen  B.  Fig.3  is  the  curve 
of  fractal  dimension  versus  SEM  magnification,  which  verifies  this  analysis.  Fig.3  shows 
that  the  fractography  cannot  be  described  by  a single  but  rather  a series  of  fractal 
dimensions.  It  also  reveals  that  the  fractography  is  not  ideally  self-affine  at  all  scales  of 
SEM.  The  transgranular  cleavage  fracture,  for  example,  shows  cleavage  facets  as 
dominant  morphology  at  low  magnification  and  river  patterns  as  dominant  morphology  at 
high  magnification,  implying  that  fractal  dimension  for  a given  mode  is  not  universal. 


Fig.3,  Fractal  dimension  Oc0r-magnification  M curves  of  samples  A and  B. 

Fig.4  shows  the  austenitic  dendrite  structures  of  HP26Cr35Ni  alloy.  Because  the 
specimens  were  taken  from  the  same  region  of  solidification,  the  distribution  and  size  of 
austenitic  dendrites  of  all  the  specimens  are  the  same  even  though  they  have  been 
subjected  to  different  thermal  shocks.  Their  impact  toughness,  however,  is  different  from 
each  other  due  to  the  formation  of  micro-cracks  during  the  thermal  shock  treatment. 
Fractal  measurements  in  Fig.4  (d)  show  that  there  is  no  significant  difference  in  the  fractal 
dimensions  of  the  specimens  under  different  thermal  shocks.  Therefore,  there  is  no  direct 
correlation  between  fractal  dimension  and  toughness  when  the  change  in  toughness  is 
caused  by  non-geometric  factors. 

The  present  experiment  has  shown  that  fractal  dimensions  of  fractography  at  SEM 
scales  are  not  a single  value  but  rather  a series  of  values  for  a given  fracture  mode.  This  is 
because  of  the  different  dominant  patterns  of  fractography  in  different  range  of 
magnifications.  It  can  be  deduced  that  when  the  magnification  is  high  enough,  the  fractal 
dimension  will  tend  to  2 because  the  micro-morphology  elements  approximately  become 
2D  surfaces.  The  irregular  geometric  structures  of  fracture  surfaces  can  reflect  the  state  of 
the  energy  absorbed  by  the  fracture  surface  due  to  the  formation  of  new  surfaces.  If  the 
change  in  toughness  is  only  result  of  the  irregularity  of  the  fracture  surface,  the  fractal  can 
be  employed. 
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Fig.4.  Microstructures,  impact  toughness  and  fractal  dimension  of  HP26Cr35Ni  alloy 
specimens  subjected  to  thermal  shock  of  different  times,  (a)  The  origin  microstructure;  (b) 
microstructure  with  10  times  of  thermal  shock;  (d)  microstructure  with  20  times  thermal 
shock;  (d)  changes  in  impact  toughness  and  fractal  dimension  with  number  of  thermal 
shock. 


5 Conclusions 


The  variation-correlation  of  fracture  surfaces  imaged  by  SEM  has  been  found  to  follow  a 
power  law  only  within  a limited  range  of  scale  of  £<£max . Beyond  this  range,  the  variation- 
correlation  tends  to  be  a constant.  The  scaling  properties  can  be  expressed  by 


Kor(e)« 


S « 8 


max 


E»S 


max 


Fractal  dimension  of  fractography  at  SEM  scales  is  not  a single  value  but  rather  a 
series  of  values  for  a given  fracture  mode.  There  exists  a positive  correlation  between 
fractal  dimension  and  toughness  of  the  materials  when  the  change  in  toughness  is  caused 
by  geometric  factors.  However,  if  the  change  in  toughness  is  caused  by  non-geometric 
factors  such  as  micro-cracks,  there  is  no  essential  relationship  between  impact  toughness 
and  fractal  dimension.  For  a given  fracture  mode,  fractal  dimension  is  not  universal. 
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